INTEGRATION OF PDES BY DIFFERENTIAL GEOMETRIC MEANS 

NAGHMANA TEHSEEN AND GEOFF PRINCE 



Abstract. We use Vessiot theory and exterior calculus to solve partial differential equa- 
tions(PDEs) of the type Uyy = F{x,y,u,Ux,Uy,Uxx,Uxy) and associated evolution equations. 
These equations are represented by the Vessiot distribution of vector fields. We develop and 
apply an algorithm to find the largest integrable sub-distributions and hence solutions of 
the PDEs. We then apply the integrating factor technique to integrate this integrable 
Vessiot sub-distribution. The method is successfully applied to a large class of linear and 
non-linear PDEs. 



1. Introduction 

The geometric study of differential equations aims to describe the local and global structure 
of the solutions as integral submanifolds of an equation manifold. This study of PDEs leads 
to the concept of Vessiot distributions, which provides the convenient formal framework to 
investigate PDEs on an appropriate jet bundle. 

In [19] the authors used exterior calculus and the Frobenius theory of foliations to show 
that integrating factors can be used to solve systems of higher order ordinary differential 
equations, and they considerably extended the class of group actions which could be thought 
of as symmetries of the system for the purpose of reduction to the quadratures. Moreover their 
methods avoided cumbersome changes of coordinates, allowing solutions and first integrals 
to be represented in the original, physical coordinates of the problem. Later this technique 
is extended to the first order PDEs in [21 [3l [U [5] . In the present paper we will extend the 
integrating factor technique to second order PDEs by using their results to integrate the 
(integrable sub-distributions of the) vector field description of the PDEs. 

The formulation of a (finite) system of differential equations as a distribution of vector fields 
on some appropriate jet space is due to Vessiot [221 1231 IM] and called a Vessiot distribution. 
The maximal integrable sub-distributions of the Vessiot distribution represent every local 
solution of the PDE system. Computing the Vessiot distribution is routine but finding the 
maximal integrable sub-distributions is not. In fact there is no known algorithm for computing 
them in the generic case. 

Solving the PDE system is equivalent to integrating a Pfaffian system or its dual vector field 
system. The Vessiot formulation of PDEs is dual to the more popular exterior differential 
system (EDS) formulation due to Cartan [6] and elaborated in [13l|T5]. In [9] the author 
pointed out that the former is in many ways computationally simpler because vector fields 
are simpler to implement than the full exterior calculus of p-forms. Recent developments in 
the Vessiot theory can be found in [lOl [201 EI]- 
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The purpose of this paper is to provide a systematic approach to solving second order Hnear 
and non-hnear PDEs in the presence of symmetries. We use a general class of symmetries 
known as solvable structures which are not necessarily of point type. To generate the solvable 
structures for each integrable sub-distributions we use the symmetry determination software 
package DIMSYM [T7] operating as a REDUCE [T2] overlay. We also use the exterior calculus 
package EXCALC [IB|. 

The paper is organised as follows. In Section 2 we introduce the key concepts of the Vessiot 
distribution. In Section 3 we give an explanation of how a reduction of order is possible 
using symmetries and present the main result of [19j which we will use to integrate the 
Probenius integrable distribution of vector fields. In Section 4 we describe, in detail, a three 
steps algorithm to find the largest integrable sub-distributions which satisfy the appropriate 
independence/projectability condition. To integrate such distributions we use the results of 
|19j . We present some examples to explain the algorithm. In Section 5 we produce an EDS 
method to solve the differential conditions in step 2 of the algorithm. We also demonstrate 
the construction of a group invariant solution of our second order PDEs. We will explain the 
methods with some concrete examples. 

2. Basic concepts and notations 

Consider a system of partial differential equations [20] of m independent variables and n 
dependent variables , 

Denote by X the space of independent variables and by U the space of dependent variables. 
The subscripts 1 < ii < . . . < < n are used to specify the partial derivatives of , where k is 
the maximum order of the system. The reader may consult [20j for more details. 

Consider the trivial bundle tt : X xU X and let / and g be two smooth sections of vr. We 
say that / and g are equivalent to order A; at x if and only if 

(9x1)91 . . . (9x™)'?™ " (9x1)91 . . . (9x™)9m 

for all m-tuples {qi ■ ■ ■ Qm) with qi + . . . + Qm < k. The equivalence class of a smooth section 
/ at a point x is the A; -jet of / at x and is denoted by j^f. The set of all A; -jets of smooth 
sections of vr constitutes the bundle of A; -jets of maps X^U and denoted by J^(X, U). The 
zeroth order jet bundle, J^{X,U), is identified with X xJJ. The A;-graph of smooth section 
/ of TT is the map j^f ■ X J^{X, U) defined by x -> j^f. Thus, the image of A;-graph of a 
section is a m-dimensional immersed submanifold of J^{X,U). 

The geometry of jet bundles is to a large extent determined by their contact structure. 
Each jet bundle J^{X,U) comes equipped with a module of differential 1-forms spanned by 

m 

2 = 1 

where / is a multi- index of order less than or equal to A: - 1 and is denoted by Q!^{X,U). 
The contact system Of'^X^U) on is also an example of an exterior differential system 
generating an ideal in the ring of differential forms on a differentiable manifold. This particular 
differential system generated by 1-forms is also called a Pfaffian system. 
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Another way to view this is that there may be a submanifold N c M where all the forms 
in the contact co-distribution are zero when pulled back to A^. That is, if z : A ^ J''{X, U) 
is any immersion whose pull-back annihilates the contact system and satisfies the transverse 
condition or independence condition i*{dx^ a ... a dx"^) 0, then i{X) is the image of some 
A: -jet. The image of the fc-jet of any function f ■ X U is an integral submanifold of the kth. 
order contact system. 

Now we can see how these constructions are related to differential equations. The study of 
second and higher order differential equations on manifolds needs to consider jet bundles on 
some underlying space. In view of the definition of a jet bundle, the A; -jet bundle J^{X^ U) 
is the appropriate framework for dealing with the systems of PDEs of a higher order of m 
dependent variables and n independent variables. We will be interested in the second order 
partial differential equations which define embedded submanifolds of J^{X,U). 

A solution of (dl) is a smooth function f : X ^ U whose fc-graph Jj defines a m-dimensional 
submanifold 5" of a smooth embedded submanifold of J^{X,U). 

3. Ideals, Symmetries and reduction of order 

We will briefly explain how the integrating factor technique in |19) allows a reduction of order 
via symmetries for higher order ordinary differential equations. We begin with some basics, 
for more details see [ill^ 171 fTHl flU] . 

Let M be some smooth manifold of dimension n. For simplicity, we suppose that all our objects 
are smooth on M. The algebra of differential forms on M, /\{M), is a graded algebra. An 
ideal I in /\{M) is an additive subgroup of f\{M) that is closed under the wedge product 
(/3* e /implies Aa e I, Ma e /\{M)). An ideal / is a differential ideal if the exterior derivative 
of every member of / is also in /. A vector field V is said to be a symmetry of an ideal if 

The kernel oi a differential form i7 is the submodule of vector fields annihilating fi, ker $7 := 
{Y e X(M) : y J i7 = 0}. A differential p-form i7 on M is simple or decomposable if it is the 
wedge product of p 1-forms, so Vl is locally simple if ker $7 is everywhere n- p dimensional, 
that is, of maximal dimension. A constraint 1-form for differential form 17 is a 1-form 
satisfying a $7 = 0, which implies Y aO -Q, MY e ker 17. An m-dimensional distribution D 
on M is an assignment of an m-dimensional subspace Dx of T^M to each point x of M. The 
corresponding co-distribution consists of the annihilators c T*M. As a submodule 
of A^(Af) the elements of are called constraint forms for D. A characterising form for 
an m-dimensional distribution D is a form 17 on M of degree (n - m) which is the exterior 
product of (n - m) constraint forms. 

A distribution D has (maximal) integral submanifolds if and only if 

(2) d0 A Q - 0, V 6 (17 is any characterising form for D) 

or 

[Y,Z]€D, yY,ZeD. 

In this case D is said to be involutive or Frobenius integrable and any characterising form 17 is 
said to be Frobenius integrable. 17 is Frobenius integrable if its kernel is Frobenius integrable 
and of maximal dimension everywhere. A function / on M is called an integrating factor for 
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the Frobenius integrable 1-form 9 if d{f9) = 0. A vector field V e X(Af) is a symmetry of 
a vector field distribution D c X{M) if CyD c D and y is a symmetry of differential form 
Q if = Afi for some smooth function A on M. If V is va. D it is said to be trivial. For 
any vector field distribution D, we say that a collection of k linearly independent vector fields 
Xi, . . . e X(M) forms a solvable structure for D if 

£x,_,(Sp{Xfc}0Z)) c Sp{Xfc}0Z?, 



£xi(Sp{X2,...,Xfc}0i^) c Sp{X2,...,Xfe}0Z^. 

DIMSYM can be used to locate solvable structures for distributions. The following results 
will be useful. 

Theorem 1. [19] Let he Frobenius integrable 1-form which is nowhere zero on some open 
subset U of a manifold M. Then fi is an integrating factor for 6 if and only if ^ - {X J 9)~^ 
for some symmetry vector field X on U with X -i9 on U. 

Proposition 2. Let 9^,9"^, 9^ e h^{M) then Sp{6'^ (9^, . . . , 6'''} is Frobenius integrable, 

A 6*^ A... A 6*^ = 0, a = l,...,A;, 
if and only if d{9^ a ... a 9^) - \ a 9^ a . . . a 9^ , for some 1-form A. 

Proof. If Sp{9^,9'^, . . . , 9^} is Frobenius integrable, then d9'^ a9^ a . . . ^9^^ - 0, and from the 
definition of exterior derivative d{9^ a . . . a9^) ^ d9^ a9'^ a . . . a9^ + . . . + (-1)^6*^ a ... a d9^ . 
It is obvious from the Frobenius conditions that {9^ a ... a 9^) divides d{9^ a ... a 9^). Hence, 
d{9^ a...a9'')^\a9^ a...a9^. 

Conversely, assume that d{9^ a ... a 9^) - X a9^ a . . . a9^ then, by using the definition of 
exterior derivative, 

d9^ a9^ a...a9'' + ... + {-1)''9^ a ... a d^^^ = a a 0^ a ... a 0^. 

By taking the wedge product with 9°", we have 

d9'' a9^ A... A9'' ^0. 

So, Sp{0^, . . . , 9^} is Frobenius integrable. □ 

We now use Proposition [2] to give a generalisation of Theorem [1] (compare with [5]). 

Theorem 3. Let Q be a k-form on a manifold M, and Sp{{Xi ...Xf^}) be a k-dimensional 
distribution on open U '= M satisfying Xi -i Q everywhere on U . Further suppose that 
Sp{{Xi . . . Xf;} u ker is Frobenius integrable for some j < k and that Xi is a symmetry of 
Sp{{Xi+i . . . Xk) u kern) for I... j. 

Put a* :- Xi J ^Xi-\ I Xf: J 0,, where Xi indicates that this argument is missing, and 

put Lo'^ — — — 7 for i - 1 . . .k, so that {u^ . . . w'^} is dual to {Xi . . . Xj^}. Then 
Xi J a* 

{i) AU)"^ A . . . Auj'' - ^ 



n{Xi,X2 . . . Xk) 
(ii) (jj^ A u}"^ A . . . A oj^ is closed. 
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Proof. To prove the first part, observe that 

We know that {Xi J Q) a {X2 J 0) = Xi J (Jl a (X2 J - {Xi J X2 J n)^. As the 2-form 9. is 
characterising form and X2 J is a hnear combination of 6 and , so that Vt a (X2 -iQ) - 
imphes that {Xi J r2) a (X2 J J]) = -(Xi J X2 J ^2)^2. Thus aco"^ ^ n{x^ Xi) ■ Clearly the 
process is inductive, and then the proof is complete. 

To prove the second part, let 



dO(X2,Xi)Aj] dn 
(j^(X2,Xi))2 ^ niX2,Xiy 
= -dn{X2,Xi) An + n{X2,Xi)dn 
= d{X2-iXi-in) An + n{X2,Xi)dn 

Using LxiCt - X2 -ida + d(X2 -i a), we have 

{n{X2,Xi)fdl-—-^ — -) = [£x^{Xi-in)-X2-id{Xi.in)]An + n{X2,Xi)dn 



n{X2,Xi) 



= -(X2 J {£x,n)An + {X2-iXi.idn)An + n{X2,Xi)dn 
= -Xi J (X2 -idn)An + {X2 J (iJ^) A {Xi J f)) + J7(X2,Xi)(io 



Now use Proposition [21 



(niX2,Xi)fdl—^^——] = -Xi j[X2J(Ka(Xi j17)a(X2 Jl^))] Al7 + 0(X2,Xi)dO 

+ [X2j(Ka(Xi j17)a(X2 a(Xi jJ^) 

= -l^(X2,Xi)dJl + 0(X2,Xi)(iJl 
= 0. 

Hence, d(cj^ a tj^) = = n{x^ Xi) j' '^'^^ proof then follows by induction. □ 

The paper by Sherring and Prince [19j extends Lie's approach to integrating a Frobenius 
integrable distribution by using solvable structures. The important results of [19] are repro- 
duced below. We illustrate these results with the straightforward case of a single second order 
ordinary differential equation. 

On some open subset of the first order jet bundle J"'^(M,M) with coordinates x,y,p = 

the second order differential equation (vector) field of the differential equation = f{x,y,p) 
is 

(3) P ^ ^ J ^ 

dx dy dp 

This has dual distribution spanned by the contact form 9 - dy - pdx and the force form 
(p - dp - dx. The 2-form = a </> on [/ is then a characterising form for ([3]). 



Proposition 4. [19] If X is a vector field whose span with T is two dimensional on U then 
the 1-form X J O is Frobenius integrable if and only if X and T are involution. 
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Corollary 5. [19! If X is a non-trivial symmetry of ([3]) then X -lO, is Frobenius integrable. 

Corollary 6. [19j If Xi and X2 are symmetries of ([3]) which commute and whose span with 
r is three dimensional then the two forms 

Xi-iQ , X2 J n 
— and — 

are closed and locally provide a complete set of two functionally independent first integrals. 

The two first integrals of this corollary, together with x, provide a coordinate chart which 
straightens out the integral curves of T. 

The next corollary to Proposition [J] shows how a solvable structure of two symmetries 
produces one closed form io^ and a second form which is closed modulo io'^, that is dio'^ = 
mod Lo^ or locally uj^ - d'y^ and - +j^d'y^. This result in [TU] is proved by using 
the Frobenius theory of foliations but we will prove it in a very simple way by using exterior 
calculus. We need Theorem O and the following lemma. 

Lemma 1. Let a,(3,df e A^(-^) be non-zero one-forms with /3 ^ Sp{df}. Then da - df a (] if 
and only ii j3 = dg mod df for some g e A°(M). Equivalently d/3^ X^df <^ P = dg mod df. 

Proof. Suppose that da - df a f5. Then it is easy to see that ■- Sp{d/,/3} is Frobenius 
integrable and so locally there exists dh such that :- Sp{d/,d/i}. Hence, (3 - m df + 
k dh, m, k € /\^ and da - df a P - df a (m df + k dh) - df Ak dh. Now = d'^a - d{df a kdh) - 
-df A dk A dh, so that dk € Sp{df,dh}. This implies A; is a composite of / and h alone. Let 
g be the partial integral of k with respect to h, then k dh - ^dh and da = df Ak dh = 
df A (^dh) - df A {^dh + ^df) - df a dg. Hence, df a (3 - df a dg which implies that 
df A (/3 - dg) - and P ^ dg mod df or P - dg + I df. The converse is trivial. □ 

Corollary 7. If Xi is a symmetry of X2 is a symmetry of Sp{Xi,T} and Xi,X2 and 
r are linearly independent everywhere then 

^ X2-iXi^n 

is closed, while 

2_ X2-in 

is closed modulo lo^. These two forms locally provide a complete set of two functionally inde- 
pendent first integrals, io^ - d-)^ and uP' - d'^'^ - X2{j'^)dj^ . By putting X2 = X2 - X2{j'^)Xi 
we then have the two commuting symmetries Xi and X2 which provide two first integrals 
and 7^ via Corollary O 

Proof. The kernel of Xi J SI is spanned by Xi and F and Xi J 17 is Frobenius integrable, since 
Sp{{Xi,T}) is closed under lie bracket. So, lo^ is closed by Theorem [TJ To see that is 
closed modulo u}'^, we consider u}^ auj"^ - Xi) using the fact d{u}^ a w^) = (see 

proof in Theorem[3|), this implies that d{uj^ a lo'^) - - d{d'y^ a uP') - dto'^ a dj^. 

If duo'^ A d'-f^ - 0, from Theorem 1 we may write duj"^ - ^ a d'y^ for some fi e so we have 
from Lemma [1] that uj"^ - dj"^ + Idj^. To find the value of /, contract this equation with Xi 
and X2, we get Xii^/"^) = 1 and / = -^2(7^). Hence, = ^7^ - ^2(7^)^7^ □ 
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Corollary [7] can be generalised to the following theorem. 

Theorem 8. [19J Let Q be a k-form on a manifold M, and let 5p({Xi, . . . , X^}) be a 
k- dimensional distribution on an open U ^ M satisfying Xj J ^fc everywhere on U. Fur- 
ther suppose that Sp{{Xj+i, . . . ,Xf^}\Jker 0,) is integrable for some j < k and that Xi is a 
symmetry of Sp{{Xi+i,.. .,Xk}\Jker n) for i^l,... ,j. 

Put fj* := Xi J ... J J ... J Xk J Q, where Xi indicates that this argument is missing and 
:= Y^-T for i - 1, . . . ,k so that {u^, . . . , u''} is dual to {Xi, . . . , Xk}- Then duj^ - 0; doj"^ - 
mod uj^] duj'^ mod a;^,Lj^; ... ]du}^ = mod uj^ , . . . , so that locally 
co^^dj\ 

uj^ ^dj^ -Xi{-f^)dj^, 

oj^ = -X2(7=')d7' - (^1(7') - X2(7')Xi(7'))'i7\ 

oj^ - dy mod dj^, . . . , ^7-'^^, 

for some 7^, . . . , 7-' e A° T*U. Also the system {o;-'"'"^, . . . ,uj^} is integrable modulo dj^, . . . , ^7-' 
and locally Q - j'^dj^ a dj"^ a ... A dj^ a uj^*^ a . .. /\uj^ for some 7° € /\^{T*U). Each 7* is 
uniquely defined up to the addition of arbitrary function 0/7^, .. . ,7*"^. 

Sherring [18] gives an algorithm for the integration of these forms which we reproduce in 
Theorem [9] below. We repeatedly use this process in our treatment of PDEs in subsequent 
sections. The integration of an exact 1-form is described in the following lemma. 

Lemma 2. Given an exact 1-form oj uji dx* on R", if u' := / (wj - Ej=if'"xi) dx^ and 
7 — E"=i then u = dj. 

Proof. Let uj = dj for some 7 as a; is exact and so Ui dx^ = ^ dx*. For z = 1, 7 = f coi dx^ + 
<p^{x'^, . . . = fj"'^ + . . . and then, 7^2 - + </'^2 = ^2 this implies that, (p^ - 

J {uj2 - 0-^2) dx'^ + . . . and 7 = cr^ + cr^ + . . . ,x"). Now assume this holds 

for i = /c, 7 = i;*Li (J* + (f)'' {x'''^^ , . . . ,x^). Then the result holds for i = A; + 1, 7 = yJI^I a* + 
. . . ,x^). The claim then follows by induction and the n*^ claim proves the theorem 
as is just the constant of integration. □ 

Theorem 9. Let oj e A(R") which is exact modulo dj^ , . . . , djK Without loss of generality 
choose coordinates x' 7', for I - and ■- x^,k - j + 1, . . . ,n. Then oj ■- 

oji dx\ put cr* := / {ui - Efc=j'+i ^Xi)"^^* foi^ ^ - j + ■ ■ ■ and 7 ■■- Er=j+i Then oj = 
mod d^y^ . . . d'yK 

Proof. Given that w is a non zero 1-form which is exact modulo dj^ , . . . , djK Clearly, by 
choosing coordinates x' = 7', / = 1, . . . , j; we then have oj is exact. Now oji - for i - j + 1, 

-f^J WJ•+ldx^■+l + (/.^■+H2;^'^^•••,2;") =aJ+i + (/)^+i(xJ+2,...,x") 

and then, 73,^+2 = o"^j+2 + '?^^ji2 - this implies that, 

<P'''- J {oj,,,-a%)dx^^' ^<lyi^\x^^\...,x-) 
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and 



Clearly, for i = j + 3, 7 = Eg+i^^'^ + 
7 - Eilj+i 0" is just the constant of integration. 




□ 



For ordinary differential equations the vector field description is one-dimensional, and hence 
integrable by dimension. For PDEs the situation is different. And so having explained the 
integration of a Frobenius integrable distribution via solvable structures we now turn to the 
issue of locating integrable sub-distributions of non-integrable distributions. 

We begin with a non-integrable distribution D of constant dimension p on an n-dimensional 
manifold M. Let be the n-p dimensional co-distribution and let D* be a complementary 
distribution with X(M) = DQD* . We aim to find the largest integrable sub-distributions of 
D, of dimension p - M say, which satisfy the appropriate independence condition. 

This is done in stages by using the following algorithm. We will try to reduce the dimension 
of and correspondingly add more 1 -forms to the distribution D. Then the integrability 
condition ^ for the augmented distribution is divided into two parts, one is algebraic and 
the other is differential 

(4) dcr" A Jlo- A Op = 0, = 1,..., n-p 



Here 0,^ — o" a ... a cr" is a characterising form for D and are 1-forms in D*^ with 
^Ip :- A ■ ■ ■ A p^^ and r^o- ^ * 0. 

Firstly, we will add a single 1-form p e D*^ to and generate the algebraic conditions 
dH) for Sp{(T^, . . . o"P,/o} to be Frobenius integrable. If the algebraic conditions can be satisfied 
then there may be integrable Vessiot sub-distributions of dimension p - 1. If the algebraic 
conditions fails then we will add two 1-forms from D*^. The process continues until we have 
an enlarged co-distribution which may be integrable, if M - p-1 then the result is trivial. The 
next step is to solve the differential conditions ([5]). This is the difficult part and in general 
the best we can do is apply an EDS process which will identify the number and functional 
dependency of the solutions rather than explicitly construct solutions. We explain this process 
in the last section. In the event that we can construct a closed form solution we can, in the 
presence of a solvable structure, integrate the Frobenius integrable sub-distribution by using 
Theorem [SI 

In the following section, we will consider the particular type of second order PDEs of one 
dependent variable and two independent variables and explain the algorithm in detail. 



In this section, we examine two types of second order PDEs of one dependent variable and 
two independent variables. The first is of the form Uyy - F{x,y,u,Ux,Uy,Uxx,Uxy) and the 
second PDEs of evolution type. 

4.1. PDEs of the type Uyy - F{x,y,u,Ux,Uy,Uxx^Uxy)- Consider a second order PDE with 
one dependent real variable u and two independent real variables x,y given by 

(6) Uyy — F i^x ,y ^u,Ux lUyiUxxiUxy) 1 



(5) 



dp" A a = 0, a = l,...,M<p. 



4. Second order PDEs 



INTEGRATION OF PDES BY DIFFERENTIAL GEOMETRIC MEANS 9 

for some smooth function F. The embedded submanifold [20] 

,y,U, Ux , Uy J Uxx ; Uxy 

is a subset of J^(IR^,M). A local solution of the PDE is a 7-dimensional locus of J^(M^,M) 
described by the map i : S J^(M^,M), i.e. 

We can study the solutions of ^ by studying the integral submanifolds N of the restricted 
contact system 

Dy ■■- Sp[6^ ■■- du - Uxdx - Uydy, 

(7) O"^ dux - Uxxdx - Uxydy, 

6^ ■■- duy - Uxydx - Fdy}, 

which project down to X c M^. If satisfies the transverse condition dx a dy\j\[ ^ and has a 
tangent space that annihilates the restricted contact structure, then i{N) c J^(M^,M) is the 
2 - graph of a solution of ([6]). So, we need to study the vector field distribution dual of the 
pulled-back contact system, that is, 

r. Q (t^ ^ 5 d d 

Dy - Sp-^ Vi-.^ — + Ux— + Uxx-;^ — + Uxy- — , 

[ ox OU OUx OUy 

d d d ^ d 

(8) V2-^ — +Uy— + Uxy^— + F—, 

oy OU OUx OUy 



dUxx OUxy 

This dual distribution is called the Vessiot distribution Dy of 



4.1.1. Algorithm for finding the largest integrable sub- distributions of Dy- 

In general the Frobenius integrability condition fails on the Vessiot distribution Dy. So, we 
will try to reduce the dimension of Dy and correspondingly add more 1 -forms from D*^ to 
the co-distribution Dy as described at the end of the Section 3. 

For this class of PDEs we need to add exactly two 1 -forms to Dy because in adding a single 
1-form (p the second algebraic condition dO^ a a (/) - implies that = 0. Furthermore, 
we cannot add more than two 1 -forms because the reduced Vessiot distribution cannot have 
dimension less than 2. 

Step 1. Add two 1-forms (1)^,(1)^ € D*^ in D^. 

We are interested in the largest integrable sub-distributions which satisfy the independence 
condition {dx Ady i= 0). So, for this purpose the coefficient of and V'^ in (p^, <p'^ should be 
nonzero. Without a loss of generality we assume that 

(9) (/.^ - aiV'^ - a2V'^ 4>^ -.^il^^ -biip^ -b2i^'^. 

Now, solve the algebraic conditions 



de"" AQeAn^^O, a = 1 ... 3. 
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Specifically, 

02 = 61, 62 = UxK + UxxK^ + aiK^^ + blK^y + UxyFuy + Fx- 

Step 2. Differential conditions: 
Solve the differential conditions ([5]) i.e. 

dcj)-^ A fig A $7^ = 0, 

d(f)^ aQq aQ^ - 0. 

Step 3. Integrate the reduced Vessiot distribution Dy^^^ by using Theorem [8l To find the 
symmetries of reduced Vessiot distribution use DIMSYM, for example. 

Remark 10. In applying the algorithm, there will exist situations when it may be difficult 
to solve the differential conditions (step 2). For such situations we will look on the alternative 
way to solve them. In Section 5, we will discuss the solution of these differential conditions 
in detail. 



We illustrate this algorithm with an example: 
Example 11. Consider the Laplace equation 



'^xx ^yy ~ 



The Vessiot distribution and co-distribution is given by equations ([7| and ^ with F - -u^a 
Now apply the algorithm step by step. By adding the two 1 -forms ([9|) in the Dy and afte 
solving the algebraic conditions dH we obtain 02 = 61, ai - -62- 

The corresponding differential constraints are 

,,1 dh db2 dbi db2 dbi db2 dbi db2 

d(p A ilg Ail^ = n^.— Uy— UxxT{ UxyT, U^yT, + U^xT, 

OU OU OX oy OUx OUx OUy OUy 

dbi db2 , dbi db2 
+ 02- bi- 61- b2- , 

OUxx OUxx OUxy OUxy 

dbi db2 db2 dbi dbi db2 dbi db2 



i2 



A ^61 A = Uy— Ux— — + — — + UxyT, Uxxl^ U, 



du du dx dy dux dux duy du 

dbi , db2 , dbi , db 
+ bi- + b2- 62- bi 



y 



^^tirj-^ d'^^^^ d'^-^xy ^'^xy 

In order to show some solutions we will choose an arbitrary value for these coordinates. For 
a particular choice of parameters we have the following solution. 

(j)^ := -ip'^^ (jp' := ^4 _ ^1 gatisfy the differential conditions and the transverse condition. 
The corresponding maximal dimension integrable sub-distribution -D^ed is generated by 

- d d d d d 

Vl - ^ + Ux — + Uxx^ + Uxi,^ + 



^ ■ -X ■ -xx ^ ^ "'xy ^ ■ n 

dx du dUx ' dUy duxy 
- d d d d d 



dy du dux duy dux 
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For the integration of reduced Vessiot distribution use Theorem [8] and the following symme- 
tries of -Drcd 

a d d d d 

— , ^2 - -X— - - — , A3 -y 



du du dux du duy ' 

d d d d .2 2sd„d„d^d 



Xi:^-xy—-y- X- , X5:=-(y -x )—-2y- 2x- 2 

OU OUx OUy OUxy OU OUy OUx ou, 

and after integration we obtain the following invariant functions 

-.^u- xux + xyuxy - yuy + ]^{uxxX^ - UxxV^ - x^y) + ^y^, 

f2 . „, „„. , ^3 . „„ .^f„2 2s 



f -Ux-yuxy-xuxx + xy, f xuxy - yuxx - Uy + -{y -x), 

1 

2' 



f'^:=Uxy-x, f-^{uxx-y)- 



The lifted solution on 5 is a common level set of these functions: 
This projects to u = ci + C2X - c^y + C4xy + \c^{x'^ - y^) + \x'^y - 

Example 12. Consider the nonlinear second order partial differential equation (generalised 
equation of steady transonic gas flow) [1^ 



Uy 

Uyy + A — + BuxUxx = 0, y * 0, 

y 

where A and B are arbitrary constants. For the sake of simplicity we choose A-1 and B - 1. 
Now apply the algorithm step by step. By adding the two 1 -forms ([9]) in the Dy and the 
solution of the algebraic conditions dH implies 02 = 61, 62 = aiUx + u\x + 

For a particular choice of parameters we find that 

y 

satisfy the differential conditions (we suppress domain considerations for clarity). The corre- 
sponding maximal dimension integrable sub-distribution D^ed is generated by 

- _ 9 d d d 



OX OU OUx ou> 

d d d , Uy, d , Uxv, d 



y 



V2 - — +^^xyT7— + [UxUxx + —JTT- - {'^xx + — )■ 



dy du dux y duy y duxy 

We now proceed to apply Theorem [8] to integrate the reduced Vessiot distribution. The 
symmetries of -Dj-ed are 

d d d ^ d d d d d 

Ai — , A2 2x— +y— +4n— + 2u^- — + ^Uy——+Uxy-z , X3 := — , 

OX ox oy on oUx oUy oUxy ou 

did d d d d d d 
X4 := logy— + -- — , X5 2u— -y— + 2ux- — + Sn^- — + ^Uxy- + 2uxx- 



du y dUy du dy dUx dUy dUxy dUx 
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After integration, we obtain the following conserved quantities: 

f - -p— (4toxx - 4ua; - y'^ul^ + {y^ul^ + 2yu^y){logu^^ + 2 logy)), 

f - log y + log Uxx - log{2uxy + yulj, 

f - i'i-QyUyUxx - l6yUxUr,y - Sy^ul^ - Sy^ul^Uxy + 4 logy), 

64m ^^^^^ ^y^'^tx + '^^y^ul^u^y - (logy)^(8/M^3. + 32y^ul^Uxy + 32y'^uly - 16) 
- logy(64yu^u^j^ + I2y'^ul^ + 32y^ul^Uxy - IGyUyUxx)), 

f - ]^\ogUxx- 

As before we can express the solution in the original coordinates. 

4.2. Evolution equations. In this section we present an additive separable solution of the 
evolution equations. Consider the following second order evolution equation 

(10) ^ F{x,y,u,Uy,Uyy), 

with an extra condition Uxy = 0. This extra condition means the solution is additively sep- 
arable. In [8J the authors say that a function has an additively separable form n(x,y) = 
v{x) + w{y) if and only it satisfies the condition Uxy - 0. 

The restricted contact distribution on is 

Dy ■■- Sp{^9^ -■- du - Fdx - Uydy, 

(P' ■- dUy - UyydyY 

The corresponding dual distribution is 

Dv Sp^ Vi:^ — + F — , V2-^ — + Uy— + Uyy — , V3 — 

[ ox OU Oy OU OUy <^Uyy 

In this case D*^ Sp{^/^"^ - dx,^'^ - dy^ip^ - duyy^ . 

Now follow the steps of the above algorithm: 
If we add a 1-form (p e D*^ to Dy and generate the algebraic conditions (jH) for Sp{^^, 0^, (/>} 
to be Frobenius integrable we find a non zero solution. So, there is an integrable Vessiot 
sub-distribution of dimension 2. The solution is 

oi = 0, 

(.-.^ _ +UyFu+ UyyFuy 

Uyy 

For the evolution equations we have found the solution of algebraic conditions by adding one 
(j) := ip^ - aiip^ - a2ip'^ in the Dy. So the process is very simple in this case we need to go to 
the step 3 directly which is a differential condition. At this stage we need an expression for 
F. For this purpose we consider the potential Burgers equation. 
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Example 13. Consider the potential Burgers equation 

_ 2 

with Uxy - 0. From ([TT]) . we have ai = 0, a2 - -2uy. This implies 

which satisfy the differential condition and transverse condition. 
The reduced Vessiot distribution is spanned by 

Dred - Sp{yi := Vi, V2 V2 ' 2UyUyyVz} . 

The final step is the integration of reduced Vessiot distribution using Theorem |8l 

The symmetries of -Dred are 

V ^ V ^ V n ^ ^ d n d 
^i-TT^ ^3-2x— + y—-Uy- 2uy 



9y ' 3u ' dx dy ^ duy duyy ' 

and after integration, we obtain 

- y- ^{ul + tiyy)"2arctanh((uy + Uyy^^Uy), 

f-U- + Uyy)x + ^ ln{Uyy) - ^ ln(u^ + Uyy), 

f - ^Hul+Uyy). 

The solution is n = C2 + xe^'^'^ - ^ In (1 - tanh^ (ye^^ - cie'^^)). 

Remark 14. It is worth pointing out the utility of our geometric approach in the investigation 
of some natural questions concerning conserved quantities and symmetries. Having solved the 
problem of finding all the maximal integrable sub-distributions of the Vessiot distribution we 
can now ask, for example, whether there exists a solution of the PDE admitting a conserved 
quantity of energy type, or whether there exists a solution simultaneously admitting certain 
specified conserved quantities. We can also ask whether there exist solutions admitting a 
particular symmetry algebra, solutions which have certain topological properties and so on. 

As an example, consider the Gaussian curvature of the graph of a solution of the PDE 

"^yy ~ -^('^) Vi^i ) Uy , Uxx ) '^xy ) • 

^ _ UxxF - uiy 



(1 + ^2+^2)2 

and suppose that this solution is the projection of a common level set of five functions f"" of 
the type we have discussed earlier. Then K is clearly constant on the graph of the solution 
if and only if dK a df^ a ... a df^ - 0, expressing that is a composite of . . . , But less 
obvious and much more useful is that K is constant on the graph if and only if dK AQgAil.^ - 0. 
Moreover, if K is not constant on the graph then the one dimensional integral manifolds of 
dK A rig A are the curves on the graph along which K is constant and their local existence 
is guaranteed. 

Further calculus with K on the solution submanifold of j2(]^2^]^^ -g gggg^g^j \)y using the 
tangent vector fields Vi, Vz- Since Vi, V2 project to the tangent directions to the graph of 
the solution, this calculus also projects. The point being that calculus on J2(M2^]k) ^jth 
functions like K is easier than on the base graph space. 
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5. Solving the differential conditions 



In this section we discuss the solution of differential conditions ([5]). We present two methods. 
The first is an Exterior Differential System (EDS) method and the second builds a solution 
invariant under a solvable structure. 

5.1. EDS method. 

From Proposition [21 we know that solving the Frobenius conditions 



for and some 1-form A. For later reference we note that, without loss, A has no 9"" or 0" 
components, and that dXAflgAQ^ - 0. This last condition means that dX only has components 
with O"' or (j)" as a factor. 

Our approach here will be to create a differential ideal of 5-forms and then, through a 
closure condition, find those ^1$ a which satisfy (|12p . 

The standard EDS reference are the books |131I15| and more details of EDS method by means 
of the inverse problem can be found in the work of [T]. We will give a brief summary of the 
method in this context. 

The EDS process for finding the 5-forms satisfying (|12p involves three steps. For the case 
of second order partial differential equations, we start with a module of 5-forms, S, and find 
the submodule of Y,^'^'^^ that forms a differential ideal. We will find (or not) our Frobenius 
integrable 5-form in this ideal. The second step is to reformulate the Frobenius condition 
in terms of an equivalent linear Pfaffian system of one forms, and final step is to use the 
Cartan-Kahler theorem to determine the generality of the solutions of the problem. 

The actual process for finding the submodule Y^^'^^^ c S that generates the differential ideal 
is the following terminating recursive procedure. Starting with the submodule T,^ ■- S find 
the submodule £ J]^ such that for all Q, e S^, dO, e (S'^), the ideal generated by S". 

We check if = S° and so is already a differential ideal. If not, we iterate the process, 
finding the submodule c c T,^ and so on until at some step, a differential ideal is found 
or the empty set is reached. If at any point during the process it is not possible to find a 
submodule containing simple forms, the problem has no solution. 

Having found 2^-', the next step in the EDS process is to express the problem of finding 
the Frobenius integrable 5-form in T.^"""^ as a Pfaffian system. Let the submodule S , be 
spanned by the set {Q^ : k - 1, . . . ,d}, and calculate 



dO" A^le ^n^^O, a = l,...,3 
d(j)" A^e aQ^^O, q = 1,...,M<2 



for 0,(f) :- (j)^ A ■■■ A (j) is equivalent to solving 



(12) 



d{i}g A Q^) - \ AQg A 



dn^ ^Ch^n- 

to find the 's which are now fixed one forms. 



.h 
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Let O = Tk^^ 1 now we are looking for those fi's in Yf^'^^^ such that dO. = A a fi, this gives: 

<^ dvk An'' + VkdO!' - TfeA A 17'' = 0, 

(13) ^{drk + rhil-rkX)An!' ^Q. 

So this leads us to the task of finding all the possible d-tuples of one forms {pk) - {pi,...,pd) 
such that 

Pfc A = 0. 

Suppose the solutions are 
with 

Pk^^'-O, 

that is, an e- dimensional module of d-tuples of 1-forms. Then equation (fT3]l can be expressed 
as 

{drk + Thil - ^kX) = -PAPk, 

(14) ^ dvk + ThCk - ^fcA + pApt = 0, 
where pA are arbitrary functions. 

At this stage, firstly the problem becomes that of solving ()14p for in terms of the unknown 
functions pA, and secondly, identifying the restrictions on the choice of these pa's. The general 
method for finding the solution for this problem in EDS is to define an extended manifold 

■- S <Si^^ with coordinates {x,y,u,Ux,Uy,Uxx,Uxy,rk,PA},k e {1, ...,d},A e {1, ...,e} 
and look for 7-dimensional manifolds that are sections over S and on which the one forms 

(15) at dvk + Vh^^ - TfcA + pApt 

are zero. To find these manifolds, ak are considered constraint forms for a distribu- 
tion on whose integral submanifolds we want. We choose a basis of forms on N, 
{ua ■= dx , dy , du, dux ■,duy , duxxi duxyiOkiT^ a} where are pulled-back basis for S, tta '= dpA 
and (Tk are defined above completes the basis. The condition that we want sections over S is 
that ai A a2 ^ 0. The Frobenius integrability condition for Du - Sp{ak} is dak = mod a^. 
But, 

(16) dak = vr^ A aj + t\^ai a 02 mod Sp{(7fe} 

where the oi a Q2 part is the torsion and vr^ are some linear combination of dp a- As dak 
expands with no dp a a dps terms, the system is quasi- linear. 

Because we want ai a 02 * on the integral manifolds, we need to absorb all the torsion 
terms into the tta - T^A-^^Oij. If any of the torsion terms cannot be absorbed, then asking for 
dak = mod Sp{(Tfc} is incompatible with the independence condition and therefore there 
is no solution. 

Once the torsion terms have been removed, so that ()16p becomes 

(17) dak^TtlAOi, 

the next step is to create the tableau fl from (|17p . from which the Cartan characters can be 
calculated allowing us to apply the Cartan test for involution. 
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In order to satisfy the integrability condition we will force all the tt^ terms to be zero. As 
a result we will obtain an over-determined system of linear algebraic equations for r^. 

We illustrate this idea with an example: 

Example 15. Consider the second order Laplace equation Uyy - -Uxx- In order to produce 
an initial ideal we solve the algebraic equation dH) giving 

02 = 61, ai- -62- 

This means that Jl^ is in the submodule, , of 2-forms generated by 

Following the EDS process above, the first step is to find a differential ideal. This is done in 
an iterative process as follows. 

Starting with T,^ - Spjil'^}, k - 1, ...,4, where fl'' - oj^ a Qg, find the submodule where 
for each ^ e S^, is in (S'^), the ideal generated by T,^. In this example we have = T,^. 
It means we already have a differential ideal and the process terminates. 

So, now having a differential ideal of dimension four, we can continue with the EDS process. 
The next step is to solve d{uj^ ^^e) - if f^oj^ ^^q, for the one forms . For this particular 
example we find that ah = 0. 

The next step is to find 4-tuples of one forms (pi, ...,^4) such that: 

(18) Pfc A w'^ A Jig = 0, 

Let pk - Pk+Pki with p^. = akiip^,pk - a^af^^- In this case, the 16 a^j are real constants satisfying 

023 = -034 - ^41, ai3 = -021 - 032, 031 =(222 -014, 042^-033 + 024, 

and the 12 d^a are arbitrary. So we have a 24 parameter family of solutions p^, A-1, ...,24, 
with p^ - for A - 13, . . . , 24 and p^ - for A - 1, . . . , 12. Following the procedure from the 
outline, we now extend 5 to a new manifold with coordinates: {x,y,u,Ux,Uy,Uxx,Uxy,rk,pA) 
and the problem can be written as that of finding the integrable distributions on with 
cjfc = where 

(19) ak-^ drk-VkX + PApt, 

along with independence condition dx Ady t Q on these distributions. Continuing the EDS 
process, set tta - dpA, Using this, a co-frame on N is {dx,dy^du,dux-,duy,duxxiduxy,(Jk-,T^A)- 
So the next step is to calculate dak modulo the ideal (cjfe) as follows: 

Taking the exterior derivative of ()19p and using this equation to remove resulting appearances 
of dvk , we have 

dak = (vTA - PA^) ^Pk - ^kdX + PAdpk, mod {ak} 

As we know that we have a 24 parameter family of solutions, so we will consider the precise 
expressions for A and dX for further simplification. Recall that A has no 6'^ components and 
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that d\ has only 6 Aip and ip Aip components. Let 

A := AiV'* and d\i Xijip^ + XiaO"'. 

dai = (tti - piX - ridXi) Aij?- + {1x2 - P2X - r2dX2) Aip"^ + (114 + vry - (p4 + pj)X - ridX^) a ifj^ 

24 

+ {-Ks - P3^ - ndXi) A ip"^ + K'^a-PaX) A +pAdp^] mod {o-fc}. 

A=13 

In this case there exist with dr^ e Sp{ip^} such that 

A = (2 - UxxAa - UxyX4)'ip^ + {uxxXa - UxyX^)il? + X^ij? + X^ip^, where A3, A4 are arbitrary, 
which can be seen by inspection from ()19p with = 0, ^ = 13, . . . , 24. 
One possible solution: 

ri = + uly, r2 = -Uxx, ra = u^y, = 1, 

and 

The distribution Sp{6^ ,9'^ ,9'^,(p^ ,<p'^} is Frobenius integrable. We now proceed to apply the 
step 4 of the algorithm. The reduced Vessiot distribution is 

d d d d d d 



^1 — o ^x o Uxx r^ "^xy o Uxx n "I" '^xy n ) 

ox OU OUx OUy OUxx OUxy 

- d d d d d d 

•'2 — ^ '^y T; Uxy 7^ Uxx T; "f^zi/ T| ^xx " 



Oy OU OUx OUy ' OUxx OUxy 

After integrating the reduced Vessiot distribution, we have 

:= e"^(ii^^cosy - n^-ysiny), e'"" (uxx sin y - Uxy cosy), 

J — U — XUx + XUxx ~ Uxx y^xy ~~ V^yi J — Uy — Uxy, J — Uxx ~ 

This projects to 

X C2 sin y - ci cos y 

u^c^ + Ciy- C5X + e — — . 

sm y - cos^ y 

5.2. The group invariance method. Li this section we discuss the group invariance ap- 
proach to solve the differential conditions on the 0": 

(20) d(/>" A fie A = 

We apply steps 2 and 3 of the algorithm together. 

Firstly we explain step 3 in this context. In applying Theorem [HI we need five symmetries 
in a solvable structure. Let Q. = VLq a be a 5-form, where i}^ comes from the algebraic 
conditions. We wish Q to be Frobenius integrable and ker Q - Sp{Fi, V2} ~ the reduced 
Vessiot distribution. Now suppose there exists a solvable structure of 5 linearly independent 
vector fields Xi, . . . ,X^. We will impose the condition Vl{Xi,X2.,X^,X4,X^) i= and the 
conditions that Xi is a symmetry of Ci, X2 is a symmetry of Xi J fl and so on down to X5 
being a symmetry of X4 j X^ j X2 J Xi J Cl. 
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These symmetry conditions are 

(21) £x,{X^^n)^X2{Xi^n), 

£xsiX^ J X3 J X2 J Xi J 17) = A5(X4 J X3 J X2 J Xi J n), 

for some smooth functions Ai, . . . , A5. Note that at this stage we have not assumed the Frobe- 
nius integrabihty of 0,. 

The way in which the symmetry conditions act as a catalyst to the solution of the differential 
conditions on can be seen more comprehensively by recalling that the sequence 

(22) n, Xi-iQ, X2 J Xi J 17, . . . , X4 J . . . Xi J 

are all simple forms and that they are all Frobenius integrable if is Frobenius integrable. 
Hence the Frobenius integrabihty of each form in the sequence is a necessary condition for 
the Frobenius integrabihty of Cl. And since each form in the sequence has one- form factors in 
Sp{9'^, 0"} their Frobenius integrabihty represents a successive simplification of the Frobenius 
integrabihty of 0. 

We will now assume further that Xi,X2,X3 are linearly independent symmetries of fig 
satisfying Xi jQg i= 0, that is, symmetries of our PDE. And we insist that the first three 
symmetry conditions of ()2ip hold. For example, this means that 

(23) Lx.n^^t^in^ (modr) 

which imposes conditions on the 4'°'. This is the essence of the technique. 

The remaining two symmetries X4,X^ in (j2ip of necessity have non-zero components, 
XY,X^ , in the kernel of Jig. If ker 17® Sp{Xi,X2,X3,X4} is Frobenius integrable then the 
one-form 

^ X4J(X3 JX2 JXi j(17gAl7^)) 

^ ^ (17eAl7^)(Xi,...,X5) 

is closed by virtue of (j2ip and it is a constraint form for 17. This closure is a necessary 
condition for the integrabihty of ker l7. We can simplify this differential condition by writing 
to as 

^e{Xi,X2,X^)Xj J 17(^ 



UJ 



ng{x,,X2,X3)n^{xY,xY) 



(It should be noted that this does not imply that .^^ is closed mod uj). 

In forcing the symmetry conditions (j2ip and the closure of uj we have found a closed form in 
17^ and reduced the two differential conditions (f^Uj) to a single condition 

(25) (i</> A 17e A 17^ = 0. 

(Note that 17^ - (p^ a cj)'^ .) This result can clearly be generalised to any dimension but in 
general it isn't true that Sp{(j)"} contains closed forms, even in this case. This is because 
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generically we don't have a solvable structure. The least we can expect is given in the following 
proposition 

Proposition 16. Suppose that satisfies 

de"" Afl^^O, a^l,...,N, 
d(j)'^ A Qg aQ^ - 0, a-l,...,m. 

Then there exist (f)'^ - I'^cf)^ , |/^| ^ and /i" such that 

(26) ^(f)'' + h29'', a = l,...,m. 

Proof. Suppose that Spj^'^jc/)"} is Frobenius integrable, 

de"" AQg An^^O, a^l,...,N, 
d(p'^ A Qg aQ^ - 0, a-l,...,m. 

This implies that 

c^^eSp{0^r}, with df A...Adf^^O. 
So, df^ - h^O"" + 1^4''^ 1 for some /i^, By relabelling the and the 0" we have 

df^ = /i^^'^ + with t 0. (since 0" e Sp{d/^}) 

Then put = /;^(/.^, a = 1, m. So that d/" = /i°6l" + (^", a = l,...,m. □ 

We close this section with the following example. 
Example 17. Consider the Gibbons-Tsarev equation 
(27) 

We follow these authors by considering /3 - 1, fi - and v - 1. Now apply the step 1 of the 
algorithm. The solution of the algebraic equation Q implies 

0,2 - bi, 1)2 - -aiUy + biUx- 

By substituting these values in the equation we have 

At this stage we have two unknown functions oi and bi . We impose the symmetry conditions 
()2ip on Cl. The symmetries of the original PDE (j27p are 

X - ^ X - ^ X - ^ 

ax aw ay 

d d . d d d d d 



d 


d 


OUy 


OUx 


d 


d 


UxX 

OUy 


dUxy 



Xi - -2l/— - 3x— - 4u— - 2Uy- UxT, + 2lixx7^ + Uxy-r, 

Oy ox OU OUy OUx OUxx OUxy 

d d d 
■- X— 2x— 2— Ux 

Oy OU OUx 
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The first three symmetry conditions ()2ip implies 

^.0, ^.0, ^.0, fl.o, ?l.o, 

ox oy ou ox ay ou 

This means that oi, 6i are independent of x, y, u. The 4*^^ and S*'* symmetry conditions 
generate further 8 conditions. We have 8 first order partial differential equations for two 
unknown functions of 4 variables. The solution of this overdetermined system is 

So, (f?- - ip"^ , (j)'^ = ip^ satisfy the differential conditions and the corresponding distribution is 
Probenius integrable. For this case we do not need to impose the further necessary condition 
on VLfj,. 

The corresponding reduced Vessiot distribution is 

- d d d d 

Vl'-^ — + UxT- + Uxx- + Uxy- , 

ox ou OUx OUy 

- 9 d 9 d 
y2-= -;:- + Uy— + Uxy- — + [UxUxy - UyUxx + — • 

oy ou OUx OUy 

After integration, we obtain 



/ -.^x + y— 

"Uxx '^xx 



._ ^2 ^ 2 2 2 2 r, : 

J — XUx — XyUxy ~ ~X Uxx + K^y^xlJ-xyUxx ~ V "^xy^xx " ^V^xy^ 

+ 2UxUxyU^^ - 2u^y - 2UyU^^ + 2Uxx): 

2 

„3 1 1 / Uxy Uxy 1 \ (.4 1 1 ^xy 

f -y+ log(n^^-^ + Uy), f -.^--loguxx, / ■ 

'^xx Uxx ^xx ^xx ^ ^xx 



This projects to 



u^-C2 + xyc4C5 - XC1C4 - yciC4C5 - ycl + c^^{y + e'^^'^'' ^'''') + ^C4(x^ + y^Cg). 



6. Summary 



In this paper, we have presented an algorithmic approach to the solution of second order 
partial differential equations of one dependent variable and two independent variables. We 
have successfully applied this algorithm to a wide class of linear and non-linear PDEs. Vessiot 
theory is used to generate the integrable sub-distributions which satisfy the appropriate in- 
dependence condition. Then a solvable structure of each integrable sub-distribution is found 
(by using DIMSYM for example). For those solvable structures of sufficient dimension, we 
integrate the reduced Vessiot distribution giving a local solution of the original system as 
the parameterised integral submanifold of the sub-distribution. These solutions appear in 
the original coordinates of the problem. Different sub-distributions correspond to distinct 
solutions. In this process we subdivide the Frobenius condition into two parts, one being 
algebraic and one differential. Solving the former is straight forward, the most challenging 
task is to satisfy the differential conditions (step 2). For the situations where the differential 
conditions are not satisfied easily, we have given two approaches to solve the problem. One is 
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an application of the standard EDS method and other involves the simultaneous use of a solv- 
able structure of symmetries to solve the conditions and to explicitly obtain the integrable 
submanifolds, giving a (generalised) group invariant solution. We give concrete examples 
throughout the paper. 

Clearly future tasks involve the generalisation to higher order and higher dimensional sys- 
tems, but we believe that the promise of symmetry reduction in the original coordinates and 
the calculus that is available for examining the shape of the solutions are exciting advances. 
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